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We discuss building models for nucleon generalized parton distributions (GPDs) H and E that 
are based on the formalism of double distributions (DDs). We found that the usual "DD+D-term" 
construction should be amended by an extra term, £E]_(x,£) built from the a/0 moment of the 
DD e(/3, a) that generates GPD E(x,£). Unlike the D-term, this function has support in the whole 
— 1 < x < 1 region, and in general does not vanish at the border points \x\ — £. 

PACS numbers: 11.10.-z,12.38.-t,13.60.Fz 



I. INTRODUCTION 

The studies of Generalized Parton Distributions 
(GPDs) [TH3] require building theoretical models for 
GPDs which satisfy several nontrivial requirements such 
as polynomiality [5], positivity [6-8 1, hermiticity [I], time 
reversal invariance [5], etc. The constraints follow from 
the most general principles of quantum field theory. Poly- 
nomiality (that may be traced back to Lorentz invari- 
ance) imposes the restriction that x n moment of a GPD 
H(x,£,;t) must be a polynomial in £ of the order not 
higher than n+1. This property is automatically obeyed 
by GPDs constructed from Double Distributions (DDs) 
[H 131 [SI IH] • (Another way to impose the polynomiality 
condition onto model GPDs is "dual parameterization" 
[TUHH]). Thus, within the DD approach, the problem of 
constructing a model for a GPD converts into a problem 
of building a model for the relevant DD. 

Double distributions F(/3, a; t) behave like usual par- 
ton distribution functions (PDFs) with respect to its vari- 
able /?, as a meson distribution amplitude (DA) with re- 
spect to a, and as a form factor with respect to the in- 
variant momentum transfer t. The factorized DD ansatz 
(FDDA) [SI [15] proposes to build a model DD F{j3, a) 
(in the simplified formal t = limit) as a product of the 
usual parton density f(f3) and a profile function h(/3, a) 
that has an a-shape of a meson DA. However, it was no- 
ticed [16] that in the case of isosinglet pion GPDs, FDDA 
does not produce the highest, (n + l) st power of £ in the 
x n moment ofH(x, £). To cure this problem, a "two-DD" 
parameterization for pion GPDs was proposed |16j . with 
the second DD G(/3, a) capable of generating, among oth- 
ers, the required £ n+1 power. It was also proposed [TB] to 
use a "DD plus D" parameterization in which the second 
DD G(f3,a) is reduced to a function D(a) of one vari- 
able, the D-term , that is solely responsible for the £™ +1 
contribution. As emphasized in Ref. [IB], one should 
also add D-term in case of nucleon distributions. The 
importance of the D-term and its physical interpretation 
were studied in further works (see Ref. |17j and references 
therein) . 

In the pion case, it was shown [18] that one can reshuf- 
fle terms between F and G functions of the F + G de- 



composition without changing the sum ("gauge invari- 
ance"). Furthermore, it was found in Ref. [19], that one 
can write a parameterization that incorporates just one 
function f(/3,a), but still produces all the required pow- 
ers up to £" +1 . A model for the pion GPD based on this 
representation was built in our paper [20] . An impor- 
tant ingredient of our construction was separation of DD 
f((3,a) in its "plus" part [f(j3,a)]+ that gives zero after 
integration over f3, and D-term part S((3)D(a)/a. For 
DDs singular in small-/? region, such a separation serves 
also as a renormalization prescription substituting a for- 
mally divergent integral over j3 by "observable" D-term. 

In the present paper, we apply the technique of Ref. 
[20] (see also [21]) for building models of nucleon GPDs 
H(x,£) and E(x,£). The paper is organized as follows. 
To make it self-contained, we start, in Sect. II, with a 
short review of the basic facts about DDs, GPDs and D- 
term, using a toy model with scalar quarks, that allows 
to illustrate essential features of GPD theory avoiding 
complications related to spin. In Sect. Ill, we describe 
the theory of pion GPD H(x,£,), presenting the results 
of Ref. [20] in a form suitable for generalization onto the 
nucleon case. In Sect. IV, we recall the basic ideas of the 
factorized DD Ansatz of Refs. [SI US]. In Sect. V, we use 
the formalism described in previous sections for building 
DD models for nucleon GPDs H(x,£) and E(x,£). 

An essential point is that two functions A and B as- 
sociated with two basic Dirac structures present in the 
twist decomposition of the nucleon matrix element do not 
coincide with H and E. In fact, A = H+E and B = —E. 
What is most important, A and B have different types 
of DD representation: A is given by the simplest (scalar- 
type) DD representation, while B is given by a more 
complicated representation coinciding with the one-DD 
parametrization of the pion case. Thus, building a model 
for H one should deal with a sum A + B, the terms of 
which have different- type DD representations. The re- 
sult of this mismatch is a term, which we call £E+(x,£) 
that is given by the "plus" part of the a/ f3 moment of 
DD e((3, a) used in parametrization for E(x, £) GPD. The 
term £DJ_(a;,£) should be included in the model for GPD 
H(x,£). However, unlike the D-term contribution, the 
function ^DjJ_(x,^) in general does not vanish both at 



the border points |x| = £ and also outside the central 
region |x| < £. 

In final section, we summarize the results of the paper. 



II. BASICS OF THEORY FOR DDS AND GPDS 

A. Matrix elements and DDs 

Parton distributions provide a convenient way to 
parametrize matrix elements of local operators that ac- 
cumulate information about hadronic structure. Vari- 
ous types of distributions differ by the nature of the 
matrix elements involved. In particular, to define 
GPDs, one starts with non-forward matrix elements 
(P + r/2\ . .. \P — r/2), with P being the average of the 
initial and final hadron momenta, and r being their dif- 
ference. In scalar case (which illustrates many essential 
features without irrelevant complications) we have 

(P + r/2|^(0){9 Ml • • • dnM0)\P - r/2) 

n-l 

= / , AnlX-Pni ■ ■ ■ "ju„_i^/i n _i + i ■ • ■ r fJ, n J 
1=0 

+ A nn {r fll . ..r Mn } . (1) 

The notation {. . .} indicates the symmetric-traceless part 
of the enclosed tensor. Since two vectors are involved, 
we have n + 1 distinct tensor structures differing in the 
number / of r factors involved. In the forward r — limit, 
only the A no coefficients are visible. Another extreme 



case is I — n, corresponding to the tensor {r Ml 



,} 



built solely from the r momentum. 

The forward r = limit corresponds to matrix ele- 
ments defining usual parton distributions f(x) as a func- 
tion whose moments produce A n $: 



f(x)x n dx = A % 



(2) 



The parton interpretation of f(x) is that it describes a 
parton with momentum xP. This definition of f{x) may 
be rewritten in terms of matrix elements of operators on 
the light cone: 



(P\^(-z/2)^z/2)\P) 



f{x)e' lx{Pz) dx + 0{z 2 ) . 



(3) 



In a general non-forward case, the parton carries the 
fractions of both P and r momenta. Note, that in the 

momentum representation, the derivative <9 M converts 
into the average fc M = (k^ + k' )/2 of the initial k and final 
k' quark momenta. After integration over k, (k) n should 
produce the P and r factors in the r.h.s. of Eq. (fTl). In 
this sense, one may treat (k) n as (j3P-\-ar/2) n and define 



the double distribution (DD) [TJ H [HI M 



(n-l)U\2 l J n 



F(/3, a)/3 n - l a L dp da = A nl (4) 



as a function whose f3 n ~ l a l moments are proportional to 
the coefficients A n i . It can be shown [TJ [31 Q~5] that the 
support region Q is given by the rhombus \a\ + \/3\ < 1. 
These definitions result in the "DD parameterization" 

(P - r/2\ij{-z/2)^(z/2)\P + r/2) 

= f F{fi, a) e -WP')-Mr*)/2 d(i da + ( 2) _ (5) 
Jn 

of the matrix clement. 



B. Introducing GPDs and D-term 

Another parametrization of the non- forward matrix el- 
ement is in terms of generalized parton distributions. In 
scalar case GPDs are defined by 

(P - r/2\tp(-z/2)tp(z/2)\P + r/2) 



e -™{P*)H{x,C)dx + 0{z 2 ) 



(6) 



and relation between GPD and DD functions is given by 

H(x, 0=[ F W, ") 6{x-0- Ca) d(3 da . (7) 
Jn 

The skewness parameter £ in this definition corresponds 
to the ratio (rz)/2(Pz). 

In the forward limit £ = 0, GPD H(x,£) converts into 
the usual parton distribution f(x). Using DDs, we may 
write 



r-l-N 
f(x) = / F(x, a) da 

-l + \x\ 



(8) 



Thus, the forward distributions f(x) are obtained by in- 
tegrating DDs over vertical lines /3 — x in the (/3, a) 
plane. As discussed above, f(x) is defined through the 
coefficients A n o corresponding to tensors without r fac- 
tors. Similarly, one can treat the A nn coefficients, corre- 
sponding to tensors without P factors, as the moments 
of another function D(a) 



D(a) (a/2) n da = A nn , 



(9) 



the D-term [TH]. From the definition of DD Q, it follows 
that 



r-l-H 

D{a) = I F(/3,a)d(3 

-i+M 



(10) 



i.e., D-term D(a) is obtained from DD F(/3,a) by inte- 
gration over horizontal lines in the {/3,a} plane. In this 



sense, one can think of "vertical" projection of DD that 
produces the forward distribution /(/3), and "horizontal" 
projection that produces D-term D(a). 
Taking the x n moment of GPD H(x,£) 

f H(x,Ox n dx = J2A nl (20 l , (11) 

J - 1 1=0 

we see that the coefficients A nn are responsible for the 
highest power of skewness £ in this expansion. 



C. DD plus D parametrization 

Parameterizing the matrix element (IT]), one may wish 
to separate the A nn terms that are accompanied by ten- 
sors built from the momentum transfer vector r only, 
and, thus, are invisible in the forward r = limit, i.e., to 
separate the D-term contribution. This can be made by 
simply using 

e -i?(Pz) = [e -i0(Pz) _ 1] + 1 (12) 

which converts the DD-parameterization into a 
"DD + plus D" parameterization 

(P - r/2\i/>(-z/2)il>(z/2)\P + r/2) 

= / [F(fi, a)}+ e -WP*)-i°(r*)/2 d p da 



(13) 



where 



+ / D(a) e~ la ( rz ^ 2 da + <D{z 2 ) , 



l-H 



[F(/3,a)] + = F{p,a)-6{/3) \ F{ 1 ,a)d 1 (14) 

-1+H 



is the DD with subtracted D-term given by Eq.(10l 
Then 



F{p,a) = [F(/3,a)]+ + 8(P)D(a) 



and 



where 



H(x,0 = H+(x,0 






(15) 



(16) 



H + (x, = / [F{j3, a)}+ S(x-/3- £a) dp da 
Ja 

= f F(P,a)\8(x-P-Ga) 
Jn L 

-6(x-£a) dp da (17) 
is the "plus" part of GPD H(x,£). 



A straightforward observation is that the x n moment 
of H+(x, £) does not contain the highest, namely the n th 
power of £, since the relevant integral 



a n [F(0,a)] + dp da 



(18) 



vanishes because the integrand is a "plus" distribution 
with respect to /3. 

For n = 0, the highest power is £°, and since the 
n = moment of H + (x, £) should not contain this highest 
power, it contains no powers of £ at all, i.e. it vanishes: 



H + {x,C)dx= / [F(P,a)} + dpda = 0. (19) 



Thus, H + (x, £) has the same property with respect to 
integration over x as a "plus" distribution 

[h(x)}+ = h(x) - S(x) j h{z)dz . (20) 

However, H + (x,£,) may be a pretty smooth function, 
without any S(x) terms. It should just possess regions 
of positive and negative values of H+(x,£) averaging to 
zero after ^-integration. 



D. D-term as a separate entity 

In the simple model with scalar quarks, one may just 
use the original DD F(/3, a) without splitting it into 
the "plus" part and the D-term. One may imagine 
that the DD F(j3, a) is some smooth function on the 
rhombus, with nothing spectacular happening on the 
P = line. In such a case, one may, of course, write 
F(P,a) = [F(P,a)}+ + 5{P)D{a), with the D-term ac- 
companied by the 5{P) function, but this term is precisely 
canceled by the ~ 8(j3) term contained in [F(/3, a)] + . 

However, if the theory allows purely i-channel ex- 
changes, then the relevant diagrams generate ~ <5(/3) 
terms not necessarily connected to other contributions. 
E.g., our scalar quarks may have a quartic interaction, 
and the f-channel loop would generate a 5(/3)tp(a) type 
contribution into F(j3, a). 

Furthermore, D term is formally given by the integral 
of F{P 1 a). An implicit assumption is that this integral 
converges, which is the case if F(/3, a) is not too singular. 
Note, however, that the integral of F(j3, a) over a gives 
f(J3), a usual parton distribution which are known to 
have a singular ~ p~ a behavior for small /3. This means 
that the /3-profile of DD F(j3, a) should be similar to 
that of f(P), and also be singular in the P — > region, 
F(P,a) ~ P~ a . The integral over P converges if a < 1. 
However, as we will see in Sec. Ill B, one may need 
the integrals involving F(/3, a)/P which diverge for any 
positive a. The integral for [F(J3,a)]+ still converges for 
a < 1, and the role of the D term in this case is to 



substitute the divergent integral 



l-|a| 



1 + M 



F(/3, a) d/3 



(21) 



by a finite function D(a) whose a n moments then give 
finite coefficients A nn . In this case, the "DD + plus D" 
separation serves as a renormalization prescription defin- 
ing the moments of DD. 

An attempt to consistently "implant" the Regge be- 
havior into a quantum field theory construction was made 
in Ref. [22], where a dispersion relation was used for an 
amplitude that has s a Regge behavior at large energies. 
For any positive a, such a relation requires a subtraction, 
which (as shown in Refs.|201 122)) results in a S((3)(p(a) 
term contributing to D{a). 



III. PION DDS AND GPDS 



A. Two-DD representation 

In fact, D-term was introduced first |16| in the context 
of pion GPDs, with pion made of spinor quarks. In that 
case, it is more difficult to avoid an explicit introduction 
of the D-term as an extra function. The basic reason is 
that the matrix element of the bilocal operator in pion 
case has two parts 

(P-r/2|^(-z/2) 7 ^(z/2)|F + r/2)| twist _ 2 

= 2PJ((Pz), (rz), z 2 ) + r„g{{Pz), (rz), z 2 ) . (22) 

This suggests a parametrization with two DDs corre- 
sponding to / and g functions [TB]. For the matrix ele- 
ment (22 1 multiplied by z^ (the object one obtains doing 



the leading-twist factorization for the Compton ampli- 
tude [23j ) this gives 

z^(P - r/2\${-z/2)^(z/2)\P + r/2) 



-i/3(Pz)-ia(rz)/2 



+ (rz)G(J3,a) 



2{Pz)F{f3,a) 
df3da +0(z 2 ). 



(23) 



Then GPDs are given by a DD representation 



H(x,£)= / [F(p,a) + ZG(p,a)]5(x-/3-{;a)dpda, 

(24) 

that involves two DDs: F{fi, a) and G(j3, a). The highest 
power £™ +1 for the x n moment of H(x,£) is given now 
by the G term, which one can separate 

G((3,a) = lG(f3,a)} + +6(f3)D(a) (25) 

into a "plus" part and D-term 

f l-\a\ 

D(a)= / G((3,a)df3 . (26) 

J-l+\a\ 



As a result, 

H(x,t) = F(x, + £G+ (.x, + sga(OD(x/0 , (27) 
where 



F{x,£)= / F((3,a)S(x-/3-ta)d(3da (28) 

Jn 



and 



G+(x,t)= / G(0,a) 8{ X -P-Ha) 

Jn L 

- 5(x - £a)] dj3 da . (29) 

The forward distribution f(x) in two-DD formulation is 
obtained from the DD F only: 



-1 — 1/31 

f[fi)= / F(P,a)da 

-1 + 1/31 



(30) 



Thus, D-term and f{x) are obtained from different func- 
tions, so the -D-term is indeed looking like an independent 
entity. 



B. One-DD representation 

Note that the Dirac index fj, is symmetrized in the lo- 

cal twist-two operators ip{lp, 9ui • • ■ dp,„}ip with the /ij 
indices related to the derivatives. Thus, one may ex- 
pect that it also produces the factor /3P M + ar^/2. In 
Ref. [24], it was shown that this is really the case. In 
other words, not only the exponential produces the z- 
dependence in the combination f3(Pz) + a(rz)/2, but also 
the pre-exponential terms come in the f3(Pz) + a(rz)/2 
combination. The result is a representation in which 

2(Pz)F(0,a) + (rz)G(P,a) 
= [2p{Pz) + a(rz)]f(fi,a), (31) 

that corresponds to 

P(/3, a) =/?/(/?, a) 

and 

G{fi,a) = af{fi,a). 

Thus, one deals formally with just one DD f(f3, a). The 
two-DD representation for GPDs (24) converts into 



H(x,g)= / (/3 + Za)f(l3,a)6(x-P-£a)d/3da 
Jn 

= x I f(P,a)5(x-0-£a)dPda (32) 

Jn 

in the "one-DD" formulation. 



The D-term in the one-DD case is given by 



IV. FACTORIZED DD ANSATZ 



rl-|a| 

D{a) = a\ f(0,a)dp 

-1+M 



(33) 



and one may write f(p, a) as a sum 



f(p,a) = [f(p,a)} + + 6(p) 



£>(«) 



(34) 



of its "plus" part 



i-M 



[f(p,a)]+ = f(p,a)-S(p) f(-r,a)dry (35) 

>/-l+|a| 

and D-term part 5(/3)D(a)/a. 

For the GPD H(x,Q, the "DD + + D" separation cor- 
responds to the representation 

H(x, = H+ (x, + sgn(OD(a:/0 , (36) 



where 



^±M1^ [ma) ,)(,•-.;-<,,, 



<J(x-£a) d^da . (37) 



Using /(/3, a) = F(/3, a)/ (3 we may rewrite 



H(x,0= / (l3 + Za)f(/3,a)5(x-P-Za)dpda 
Jn 



F(/3,a) 6(x - /3 - £a) dp da 



aF(P,a) 



S(x — P — £a) — S(x — £a) 



dp da 



in P 
sgn(OD(x/0 

: F DD (x,t)+SFl(x,0+sgn($D(x/0 , (38) 



where 



F DD (£,£)= / F(P,a)5(x-P-£a)dpda (39) 

is GPD constructed from DD F(P,a) by the same for- 
mula as in scalar case. Another term 

Fl(x,Z) = j^(jF{p,a)\ 5(x-P-£a)dpda (40) 



is a GPD built from the "plus" part of the DD 
aF(P,a)/p. The latter, of course, may be written as 
G(P, a), but in the spirit of the one-DD formulation, one 
may wish to express the results in terms of just one func- 
tion F(p,a). 



In the forward limit £ = 0, GPD H(x,£) converts into 
the usual parton distribution f(x). In the one- DD for- 
mulation, we may write 



f(x) =x I f(x,a)da 

-1+1*1 



(41) 



Thus, the forward distributions f(x) are obtained by in- 
tegrating over vertical lines p = x in the (P,a) plane. 
For nonzero £, GPDs are obtained from DDs through 
integrating them along the lines P = x — £a having l/£ 
slope. The reduction formula (41 ) suggests the factorized 
DD Ansatz 



f(P,a) = h(0,a)f(P)/P 



(42) 



where /(/3) is the forward distribution, while h(P,a) de- 
termines DD profile in the a direction and satisfies the 
normalization condition 



1-1/31 



h(P, a) da = 1 



(43) 



-1+1/31 



The profile function should be symmetric with respect 
toa-> —a because DDs f(P,a) are even in a [TBI |2"5] . 
For a fixed P, the function h(0, a) describes how the 
longitudinal momentum transfer r + is shared between 
the two partons. Hence, it is natural to expect that 
the shape of h(P, a) should look like a symmetric me- 
son distribution amplitude (DA) <p(a). Since DDs have 
the support restricted by \a\ < 1 — \P\, to get a more 
complete analogy with DAs, it makes sense to rescale 
a as a — (1 — |/3|)7 introducing the variable 7 with p- 
independent limits: — 1 < 7 < 1. The simplest model 
is to assume that the 7-profile is a universal function 
5(7) for all p. Possible simple choices for 5(7) may be 
J (7) (no spread in 7-direction), |(1 — j 2 ) (characteristic 
shape for asymptotic limit of nonsinglet quark distribu- 
tion amplitudes), ff(l — 7 2 ) 2 (asymptotic shape of gluon 
distribution amplitudes), etc. In the variables j3, a, these 
models can be treated as specific cases of the general 
profile function 



h { - N \p,a) = 



T(2A + 2) [{l-\p\f-a 

22JV+i r 2(JV + l) (I _ |/3|)2JV+1 



21 AT 



(44) 



whose width is governed by the parameter N . 

To give a graphical example, we show in Figjllthe sim- 
plest GPD F DD (x,g) Q built from the model 



F(P,a) = f(p)h^(p,a) 



with forward distribution 



r°«(x) = (1- xy 



(45) 



(46) 



and N = 1 profile function (analytic results for 
/ mod (a;) = (1 — x) 3 x~ a and N = 1 profile may be found 



F DD (x,i) 




FIG. 1. GPD F DD {x, f) for £ = 0.3. 



in Refs. 26-28J). The model forward function was cho- 
sen in the form reproducing the x — > 1 behavior of the 



nucleon parton distributions and the 



„-0.5 



Regge be- 



havior of valence part of quark distributions for small x, 
which was taken for simplicity, though the GPD shown 
corresponds to the C-even component (the full function 
is antisymmetric in x, and only the x > part is shown). 



In the forward limit, we have 

A(x, 0) = H(x, 0) + E(x, 0) = f(x) + e(x) (50) 

and 

B(x, 0) = -E(x, 0) = -e(x) . (51) 

These reduction formulas suggest the model representa- 
tion 



and 



a(p,a) = f(P,a) + e(P,a) 

e((3,a) 



HP, a) 



P 



(52) 



(53) 



Because of possible singularity of e(/3, a)/ (3 at /? = 0, we 
write it in the "DD + + D" representation: 

KM = -(<^\ +6W°M t (54) 



P 



where D{a) is the D-term. 



n 



B. General results for GPDs 



V. NUCLEON GPDS 

A. Definitions of DDs and GPDs 

In the nucleon case, for unpolarized target, one can 
parametrize 



W(-z/2)^(z/2)|p)| twist _ 



(47) 



-i/3(Pz)-ia(rz)/2 



+ 



n 
u(p')u(p) 
2M N 



u(p')fu(p)a(P, 



[2P(Pz) + a(rz)]b(j3,a) 



df3 da +0(z 2 



Here, the functions a, b are DDs corresponding to the 
combinations A = H + E and B = —E of usual GPDs 
H and E (see Ref. |28J. These GPDs may be expressed 
in terms of relevant DDs as 



A(x,0 



a(/3, a) S(x — j3 — £a) d/3 da 



(48) 



and 



B{x, = x I b{/3, a) 5{x-fi- fa) d/3 da . (49) 
Jo. 

Notice that we have two different types of relations 
between GPDs and DDs: A(x,£) is obtained from its 
DD a(j3, a) just like in the simplest scalar case, while 
B(x, f ) is calculated from 6(/3, a) using the formula with 
the one-DD representation structure. The difference, of 
course, is due to the factor [2j3(Pz)+a(rz)] in the 6-part. 



As a result, we have 

H{x^)=A(x^) + B{x 7 

[f(0, a) + e(/3, a)} S(x - j3 - fa) d/3 da 



(55) 



P J. a 



S(x — ft — fa) d/3 da 



= F DD (x,0 + E DD (x,0-E+{x,0+sgn(OD(x/0 
where 



Fdd(x,£)= / f(/3,a)8(x-p-Za)dj3da (56) 



and 



E DD (x,0= / e(p,a)6(x-P-£a)dPda (57) 
Jn 

are build from DDs by simplest formulas not involving a 
division by /? factors, while 



E+(x,Q f e(P,a) 
x Jn P 



6(x — P — fa) — S(x — fa) 



dp da 



e(P,a) 

P 



6{x - P - fa) dp da (58) 



has the structure of a one-DD representation. Since 
E + (x, f)/a? is built from the "plus" part of a DD it should 

satisfy 



f E + {x,i)dx= f 
J-i Jn 



e(p,a) 



dpda = 0. (59) 
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FIG. 2. GPD E+(x, £) for £ = 0.3. 



Being (for C-even combination) an even function of x, 
the function E+(x,£)/x obeys 



1 E+(x,£) 



dx = 



(60) 



C. Modeling GPDs 



To illustrate the structure of E + (x,!;), we show it in 
Figj2] using the model based on 



e(P,a) = e(P)h ( - 1 '>(P,a) 



(61) 



with N = 1 profile function and the same forward dis- 
tribution e(x) = (1 — x) 3 /y/x that was used to model 
Fdd above. Again, we have in mind the C-even, 
quark+antiquark part of the distribution, and valence- 
type functional form is used to simplify the illustration. 
One can see that E+(x, f ) is a regular function, and van- 
ishing of E + (x, £)/x integral is due to compensation over 
positive and negative parts (see Fig|3| rather than be- 
cause of subtraction of a S(x) term. 

In a more realistic modeling, one should adjust nor- 
malization of e(x) to reflect its relation to the anomalous 
magnetic moment. Also, the fits of the nucleon elas- 
tic form factors [5S] suggest for e(x) a higher power of 
(1 — x). However, our aim while showing the curves in 
the present paper is just to illustrate the qualitative fea- 
tures of various GPD models, so we will stick to the same 
generic forward function both for f{x) and e(x). 



The function E + (x,t;) may be displayed as 
e(j3,a) 



E + (x,Z)=x f 
Jn 



[6(x-P- £a) - S(x - fa)] d/3 da 



■a P 
e(P, a) 8(x — P — £a) dp da 



+ f / - e(P, a) [5(x - £ - f a) - S(x - fa)] dp da 
Jn P 

= E DD (x,£)+£ J (^e{P,a)J S(x - p - fa) dp da 

= E DD (x^)+^E 1 + (x,0, 
where 

' a 



(62) 



El(x,0= / (=e( .i.m 



5{x ~ P ~ £a) dp da 



(63) 



Since E]_(x,t;) is built from the "plus" part of a DD, 
its x-integral from —1 to 1 is equal to zero, but in fact 
it vanishes also for a simpler reason that E±(x, f) is an 
odd function of x. So, in this case, we cannot make 
any conclusions about the magnitude of the x-integral of 
E\{x,(,) from to 1. 

Summarizing, GPD E+ is obtained from the naive 
Edd function by adding to it the fi£J_(a;, f) term, which 
results in a rather nontrivial non-monotonic behavior of 
the E + function. To get the full GPD E, one should 
subtract also the D-term contribution: 

E(x, f ) =E + (x, f ) - sgn(f )£>(x/f ) (64) 

= E DD (x,t)+£El(x,Z) ~ sgn(f)D(a;/f) . 

For GPD H, we then have 

H(x,0 = F DD (x,0-^El(x,0+sga(OD(x/0- 

(65) 

Now one should subtract f£^_(;r,f) from the naive Fjjjj 
function and then add the D-term contribution. 
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FIG. 3. GPD E + (x,£)/xior£ = 0.3. 
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involves the integral 



P 



e(/3, a) I df3 da 



(67) 



that vanishes because the integrand is a "plus" distribu- 
tion with respect to j3. Hence, ££^(x,£) term also can- 
not produce the £™ +1 contribution for the x n moment of 
H(x,£). Such a term is produced by the D-term only. 



E. Comparison with "DD plus D-term" model 



FIG. 4. Function £E]_(x,£) for £ = 0.3. 
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FIG. 5. 
£ = 0.3. 



Model nucleon GPD H(x, £) (without D-term) for 



Comparing this result with the pion case for which 
H(x, = F DD (x, + £Fl(x, + sga(0D(x/0 , 



(66) 



we see that the structure of the pion GPD H + is similar 
to that of the nucleon GPD E + : the term £F}(x, £) is 
added to Fjjd(x,£,) rather than subtracted. However, in 
case of the nucleon GPD H, the extra term is built from 
the second nucleon DD e(/3, a) rather than from f(f3, a), 
and it is subtracted from Fdd(x,^) rather than added to 
it. 



The usual "DD plus D-term" model in the context 
of the present paper corresponds to "Fjjd plus D-term" 
combination, i.e. modeling nucleon GPDs without sub- 
tracting the £.&J_(:e,£) term when modeling H(x,£,), (or 
adding it when modeling E(x,t;)). 

In a sense, our new model results from the old 
"DD plus D" model by substituting sgn(£)D(x/£) with 
-ZEX(x,Z)+s&i(OD(x/£). 

Since the D-term is fitted to data, one may wonder 
if adding £_EjJ_(a;, £) may be absorbed by redefinition of 
the D-term. However, there are important qualitative 
differences between E+(x,£) and D(cc/£). First, the sup- 
port region of E\{x, £) is not restricted to the segment 
| as j < £. Furthermore, existing models of D(x/£) assume 
that it is a continuous function that vanishes not only 
outside the central \x\ < £ region, but also at the bor- 
der points |jc| = £ (otherwise, GPDs H and E would be 
discontinuous at the border points, and pQCD factoriza- 
tion formula for DVCS would make no sense). As we 
have seen, E+(x,£) is a continuous function of x in the 
whole \x\ < 1 region, and it is not vanishing at the border 
points |a:| = £. 

Thus, the most apparent difference between the two 
models is that the value of H (£, £) , the GPD at the border 
point, in the new model is different from that given by 
GPD F DD (£,g) built solely from DD f(fi,a) related to 
the usual forward parton density /(/3). Furthermore, this 
difference is determined by DD e(j3, a) that is related to 
GPD E(x,£) invisible in the forward limit. 



VI. SUMMARY 



D. Poly normality 

Taking the x n moment of H(x, £) in this construction, 
we note that the Fdd{x, £) term produces only the pow- 
ers of £ up to £™. Next observation is that the highest, 
namely n power of £ in the x 11 moment of E\{x, £) 



Summarizing, the model for GPD H proposed in this 
paper differs from the "old-fashioned" DD+D model 
by an extra — £_Ejj_(x, £) term constructed from the DD 
e((3,a) corresponding to the GPD E(x, £). The inclu- 
sion of such a term modifies the original DD-based term 
Fdd(x,^) at the border points \x\ — £ and outside the 
central |x/£| < 1 region, which may have strong phe- 
nomenological consequences. 



ACKNOWLEDGEMENTS 

I thank H. Moutarde and A. Tandogan for discussions, 
and C. Mezrag for correspondence. 

This work is supported by Jefferson Science Asso- 
ciates, LLC under U.S. DOE Contract No. DE-AC05- 
06OR23177. 



[1] D. Mueller, D. Robaschik, B. Geyer, F. M. Dittes, and 
J. Horejsi, F ortschr. Phys., 42, 101 (1994), |arXiv:hep-| 

ph/9812448] 

X.-D. Ji, |Phys. Rev. Lett., 78, 610 (1997)| |arXiv:hep 



[3] 
[4] 






[5 



[6 



T8J 
[9] 

M 



Tit 

[12 
[13 



14 



TT5 
[16 



ph/9603249 

A. V. Rad yushkin, Ph ys. Lett., B380, 417 (1996) 

|arXiv:hep- ph/96043T7p 



J. C. Collins, L. Frankfurt, and M. Strikman, Phys. Rev 



D56, 2982 (1997), arXiv:hep-ph/9011433. 

X.-D. Ji, |J. Phys., G24, 1181 (1998)| |arXiv:hep- 



(1998) arXiv:hep-ph/9711371 

B. Pire, J Soffer, and O. Teryae v, |Eur. Phys. J., C8,| 



103 (1999), arXiv:hep-ph/9804284. 



A. V. Radyushkin, Phys. Rev., D59, 014030 (1999) 
arXiv:hep-ph/9805342 



A. V. Rad yushkin, ]Phys. Lett., B385, 333 (1996) 
arXiv:hep-ph/9605431 



M. V. Polyakov and A. G. Shuvaev, (2002), arXiv:hep 



ph/0207153 



M. V. Polyakov, (2007), arXiv:071 1.1820 [hep-ph] 

M. V. Polyakov, |Phys. Lett., B659, 542~ (2008) 

arXiv:0707.2509 [hep-ph] 

K. M Semenov-Tia n-Shan sky, |Eur. Phys. J., A36, 303| 



(2008) , arXiv:0803.2 218 [hep-ph]] ~ 

M. V. Polyakov and K. M. Semenov-Tian-Shansky, Eur. 

•^rz = : — — rrr- — , . i i ^^r _ _ . . — _ _ _ . — tz r-^ ' 



Phys. J., A40, 181 (2009), arXiv:0811.2901 [hep-ph] 

A. V. Radyushkin, |Phys. Lett" B449, 81 (1999) 

|arXiv:hep- ph/9810466T~ 

M. V. Polyakov and C. Weiss, [Phys. Rev., D60, 114017| 



(1999), arXiv:hep-ph/9902451 
[17] K. Goeke, M. V. Polyakov, 



and M. Vanderhaeghen , 



Prog. P art. Nucl. 
ph/0106012 



Phys., 47, 401 (2001) arXiv:hep 



[19 



[20 



ph/9807358 [21 

A. P. Martin and M. G. Ryskin, |Phys. Rev., D57, 6692 

[22 



[23 



1 20 



O. V. T eryaev, P hys. Lett., B510, 125 (2001)||arXiv:hep 
ph/0102303 



A. V. Belitsky, D. Mueller, and A. Kirchner, Nucl. Phys., 
B629, 323 (2002)] |axXTv:hep-ph/0112108" 



A. Ra dyushkin, Phys. Rev., D83, 076006 (2011) 



arXiv:1101.2165 [hep-ph] 

A. V. Radyushkin, Int.J.Mod.Phys.Conf.Ser., 20, 251 



Londergan, 
Polon., B40, 



and F. J. 
2193 (2009), 



Nucl. Phys., B311, 541 



(2012)] 

A. P. Szczepaniak, J. T 

Llanes-Estrada, Acta Phys. 

arXiv:0707.123 9 [hep-ph]] 

I. I. Balitsky and V. M. Braun 

(1989) 

A. V. Belitsky, D. Mueller, A. Kirchner, and A. Schafer, 

Phys. Rev., D64, 110002 (2001), arXiv:hep-ph/0011314 

L. Mankiewicz, G. Piller, and T. Weigl,|Eur.Phys.J ., C5, 

119 (1998)j [arXiv:hep-ph/9711227 [hep-ph 



I. V. Musatov and A. V. Radyushkin, Phys. Rev., D61, 



" 074027 (2000), arXiv:hep-ph/9905370 
[27] A. V. Radyushkin, (2000), arXiv:hep -ph/0101225 
[28] A. V. B elitsky and A. V. R adyushkin, |Phys. Rept., 418, 

1 (2005) [ |arXiv:hep -ph/0504030 
[29] M. Guidal, M. Polyakov, A. Radyushkin, and M. Van- 



derhaeghen, Phys.Rev., D72, 054013 (2005) arXiv:hep 
ph/0 410251 [ hep-ph] | 



